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Abstract 

A review is made of the basic tools used in mathematics to define a 
calculus for pseudodifferential operators on Riemannian manifolds en¬ 
dowed with a connection: existence theorem for the function that gen¬ 
eralizes the phase; analogue of Taylor’s theorem; torsion and curvature 
terms in the symbolic calculus; the two kinds of derivative acting on 
smooth sections of the cotangent bundle of the Riemannian manifold; 
the concept of symbol as an equivalence class. Physical motivations 
and applications are then outlined, with emphasis on Green functions 
of quantum field theory and Parker’s evaluation of Hawking radiation. 


1 Introduction 

In the course of studying partial differential equations on M"", one discovers 
one can consider operators whose action is defined, at least formally, by the 
integral 

(Pn)(x) = (27r)"t f e‘^'^p(x,0h(0dG (1) 

JR" 

where d^ is the Lebesgue measure on M"’, p{x,^) is the amplitude of the 
operator P, ^ ■ x = (x,^) is its phase function, and 

u{^) = (27r)~2 f e~^^'^u{x)dx ( 2 ) 

JR" 

is the Fourier transform of u. These operators are said to be pseudodifferen¬ 
tial and form a class large enough to contain the differential operators, the 
Green operators and the singular integral operators that are used to solve 
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partial differential equations. They also contain, for each elliptic operator, 
its parametrix, i.e. an approximate inverse up to an operator of lower order 
(see below). 

To understand how pseudo-differential operators can be used in solving 
inhomogeneous partial differential equations, let us begin by considering, for 
n > 3 and a smooth function / with compact support, i.e. / € 
the inhomogeneous equation Au = /, where 

^^dxl 

k=l ^ 

is minus the Laplacian (with our convention, A is defined in the standard 
way, but the Laplacian has symbol given by = I^P)- Upon taking 

the Fourier transform of both sides, one finds 

- = /(O, (3) 


and hence, for the inverse operator Q of A, or fundamental solution, we can 
write [UE] 


u{x) = {Qf){x) 


_£(i^ f m 

4vrt L\x-y\-^y 


(4) 


If, instead of the Laplacian on M"’, we deal with a general partial differen¬ 
tial operator with constant coefficients which can be written as a polynomial 


P = p{D), where D = ..., ^ j , so that 


p{D)u = / € Co°°(]R"), 


(5) 


the solution can be formally expressed as 

u{x) = {Qf){x) = (27r)"t f (6) 

JR" 

where the amplitude q{^) is the inverse of the symbol p(^) of p{D), and the 
integration contour must avoid the zeros of p(^). However, if the operator 
P has variable coefficients on a subset U of M”, i.e. 

P = p{x,D)= ^ aa{x)D^, aa£C°°{U), (7) 

a:\a\<k 

one can no longer solve the equation Pu = / by Fourier transform. One 
can however freeze the coefficients at a point xq G U and consider P as a 
perturbation of p(xo, D), which is hence a differential operator with constant 
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and 


coefficients. In this way the amplitude of P reduces to q{^) = , 

if we let xq vary in U, we obtain the approximate solution operator 


{Qf)ix) = i2Tr) 2 [ e‘^'^g(x,0/(Od? 

JR'* 


( 8 ) 


with amplitude q{x,^) = pix,^)~^, for x G C/, / G 

However, to solve the inhomogeneous equation Pu = /, we do not strictly 
need the full inverse operator or fundamental solution, but, as we said before, 
it is enough to know a parametrix, i.e. a quasi-inverse modulo a regularizing 
operator. One can provide a hrst example of parametrix by reverting to 
the study of constant coefficient operators. A parametrix Q can then be 
constructed by choosing its amplitude [T] as 

.({) = f}. (s) 

where x(^) is a suitably chosen function which is identically zero in 

a disk about the origin and identically 1 for large In this way the integral 
formula 

(Q/)(x) = (27r)-t [ (10) 


is not affected by the convergence problems that would be met if the ampli¬ 
tude were taken to be just H Q is the integral operator in (fTOjl it is 

no longer true that PQ = QP = I, but we have [T] 


PQf = f + Rf, feC^^iU), 

(i?/)(x) = (27r)"t / r(x-0/(0dC, 
JR" 


( 11 ) 

( 12 ) 


where the Fourier transform of r is y — 1. Thus r is a smooth function and R 
turns out to be a smoothing operator. Such a class of smoothing operators 
is fully under control, and hence a parametrix Q serves just as well as a full 
fundamental solution, for which R vanishes identically. 

In the following we will see how to extend the theory of pseudodiffer¬ 
ential operators on M” to more general pseudodifferential operators defined 
on compact manifolds. In particular, we will outline some basic symbolic 
calculus for such operators. In the last two sections of this paper, some 
applications to physics will be shown. 


2 Pseudodifferential operators on manifolds 

First, note that the material in the appendix can be re-expressed by saying 
that a pseudodifferential operator A acting on functions in C“(M’^) has a 
symbol given by 




crA{xo,C) = 




(13) 
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x=xo 







and hence formal application of Taylor’s formula yields 




E 

ki .. .kn 


-h^n) 

hl...knl 


dki 

w 


d^^aA\ ( 

d^t J 



(14) 

the sum being taken over all values of the multi-index k = (ki,kn)- The 
symbol of the product AB of two pseudodifferential operators A and B is 
then defined by 


aABix,0 = e-^^-^ABe^^-^ = 


= E 

k\ ...kn 


- \~kn) 

kil...knl 


dkl 


dit J 



(15) 


Note here the crucial role played by the function 

n 

lixo,C, x) = ^ ■ (x - Xo) = ~ ^o) (16) 

i=i 

which, for each xq, is linear in x and Its derivative with respect to x is 
while its derivative with respect to ^ is x — xq. 

On going from M”’ to compact manifolds, we look for a real-valued func¬ 
tion 


I : T*M X M ^ M 

{v, x) !->• l{v, x), 

which generalizes the function ()16l) . Linearity in ^ becomes linearity in v on 
each fiber of the cotangent bundle of M, but linearity in x has no obvious 
counterpart. However, if there exists a connection V on r*M, linearity at 
xo can be defined by stating that, for all integer k > 2, the symmetrized 
fe-th covariant derivative vanishes at xq- The desired linear function is then 
a real-valued function I G C°°{T*M x M) such that the image l{v,x) is, for 
fixed X, linear in each fiber of T*M, and such that, for each v G T*M, 


d J x)\x=TY(^y'j 


V if A; = 1 
0 otherwise. 


(17) 


With this notation, is the symmetrized /c-th covariant derivative with 
respect to x, while vr is the projection map vr : T*M M. 

Once a connection V is assigned, a definition of symbol of a pseudodif¬ 
ferential operator A on C°°{M) is provided by [3] 


CTAiv) 




x=7r{v) 


(18) 
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which is a generalization of formula (|13p . However, the function l{v,x), 
whose existence will be proved in the next section, is not uniquely deter¬ 
mined by the linearity conditions above. Therefore, different functions I 
would lead to different symbol maps. On the other hand, it can be proved 
[3] that the difference between any two symbol maps corresponding to dif¬ 
ferent choices of functions I belongs to a certain class of functions, therefore 
the symbol of a pseudodifferential operator will be actually defined as an 
element of a quotient space suitably defined. 


3 Existence theorem for the function I 


Following our main source [3], we are now going to prove that there exists 
a function I E C°°{T*M x M) with the properties listed above. Indeed, we 
have the following result. 

Proposition 3.1 ([3]). There exists a funetion I E C°°{T*M x M) sueh 
that l{-,x) is, for eaeh x E M, linear on the fibers ofT*M and such that, 
for eaeh eotangent vector v E T*M, eg. holds. 

Proof. The desired function I is first constructed locally. If ?7 is a coordinate 
neighbourhood in M with local coordinates x* and an m-th order covariant 
tensor its covariant derivative is given by 

d 

V 


E 


tt 


r, 


. ..trn '> 


(19) 


where is the standard notation for Christoffel symbols. Thus by induc¬ 
tion, for any scalar function /, one has 


d^f V- . , 

dx^^---dx^k ^ 


( 20 ) 


where j are multiindices and y’s are polynomials in the derivatives of Christof¬ 
fel symbols. This implies that, for k > 1, condition (IIZI) is equivalent to each 
term 

dH 

... dx'^k 

being equal, at 'k{v), to some linear combination of lower-order derivatives. 
To sum up, starting with the requirements 


/(u,7r(u)) = 0, 


dl{v,x) [ ^ \ 

X=n(v)~''\dx0 ' 


( 21 ) 


specifies what 


x=7r{v) 


must be, in order eq. (1171) to hold. Borel’s theorem 


ensures that there exists a C°° function having partial derivatives arbitrarily 
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prescribed, and the proof shows that I may be chosen to be both linear and 
in V [3]. 

Having established that, for each coordinate neighbourhood Ui in M, 
there exists an li G C°^{T*Ui x Ui) with the desired properties, we can take 
finitely many Ui covering M with a partition of unity given by functions 
ifi € C^{Ui), and yet other smooth functions with compact support V’i € 
C^{Ui) equal to 1 on a neighbourhood of the support of The function 


l{v,x) = '^ipi{TT{v))li{v,x)ipi{x) 
i 

is then globally defined and satisfies all requirements. 


( 22 ) 


□ 


4 Analogue of Taylor’s theorem 

For a smooth function / on M, is replaced by V^/(xo). To obtain 

the analogue of x — xq note that, since for fixed xq and x G M the function 
l{v, x) is linear for v G we may think of it as an element of the tangent 
space at xq- We can instead regard Z(-,x) as a vector field on M, so that 

/(•, x)^ = Z(-, x) (g) • • • ®/(•, x) (23) 

is a symmetric k-th order contravariant tensor field, and hence 

V^f{xo)-l{xo,x)’^ 

is defined and, by virtue of symmetry of /(xq, x)^, it coincides with d^f{xo) ■ 
l{xo,x)^. A basic theorem [3] holds according to which, for each point 
Xq G M and each integer N, one has 

N 


d’^fixo) = <9^ V ^V”/(xo) • /(xo, X 
ni 


n=0 


(24) 


X=XQ 


for k < N. 


5 Torsion and curvature terms in the symbolic cal¬ 
culus 

In the symbolic calculus, one encounters frequently the unsymmetrized co¬ 
variant derivatives 

M ’ 

X=7T{V) 

which turn out to be polynomials in the torsion tensor TU and curvature 
tensor R^ij^- Indeed, by virtue of the Ricci identity, the difference of second 
covariant derivatives of an m-th order covariant tensor r,-. 


= VH{v,x) 


is given by 




E 


...In 


^ii,jk ( 26 ) 


jk- 
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Thus, for any permutation a of 1,..., fe, the difference 


is a sum of terms each of which is a product of {k — l)st or lower-order 
covariant derivatives of / and covariant derivatives of T and R followed by 
contraction. In particular, if / = /, one has at x = Tr(v) 

^ = 0- (27) 

a 

In particular, one finds, using the Einstein summation convention, 

(V l)ij = l-ij = —VpTkj, (28) 

and, in the case of vanishing torsion, 

{y^l)ijk = l,ijk = (^^) 


6 The derivatives and 


Given a function a € C°°{T*M), one defines D^a to be the k-th derivative 
of a in the direction of fibers of T*M. Thus, for 7r(u) = xq, think of a as 
a function on the cotangent space at xq, and take its A:-th derivative D^a 
evaluated at v. This is a fc-linear function on T*xo and may be identified 
with an element of the tensor product ®kTxo- This means that D^a is a 
contravariant /c-tensor, and it is the analogue of , but of course k is an 

integer in and a multi-index in 

The covariant derivatives act on C°°{M), and to define their action 
on (T G C°°{T*M) we set [S] 


VV(x) = VV(d,,Z(x,x)) 


x=7r(v) 


(30) 


Although the function l{v, •) is not unique as we said before, all its derivatives 
are determined at '/r(u), and hence (l30|) dehnes unambiguously as a 
covariant A:-tensor. The mixed derivatives may also occur and are 

dehned by 

V'^D^aiv) = u(d,/(u,x))U=^(,) . (31) 

This is a contravariant (resp. covariant) j-tensor (resp. A:-tensor), or tensor 
of type {j, k). For example, given a Riemannian manifold {M,g), if a is the 
squared norm of v, i.e. 


a{v) = \v\^ = g^'^ViVj, 


(32) 
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one has 


(33) 


{Day = 25 *%,, {D^af = 

Moreover, since 

a{dl) = (34) 

and both and g^^ have vanishing covariant derivatives at vr(r;), one finds 

Vcj = 0, (35) 


while, by virtue of ([29]), 

(VV)„ = (««, + «”/,). (36) 

and 

i^^Da)\, = y,[R«, + R\-,). (37) 

Note also that a{v) = \v\‘^ is the symbol of the Laplacian on a Riemannian 
manifold, if the connection with respect to which the symbol is taken is the 
Levi-Civita connection. 

7 The symbol as an equivalence class 

Let (M” X M™), for 1/2 < p < 1, be the space of functions a G C'°°(M"' x 
MD) such that, for all multi-indices j and k, 

= O ((1 + |C|)-^l^l+(i-^)bl) , (38) 

uniformly on compact x-sets. 

On going from M"" to a (compact) manifold M, the spaces Sy{M x 
MD) and Sy{T*M) consist of functions satisfying (I38p in terms of local 
coordinates. If it is sufficiently clear what the underlying space is, one 
writes simply Sy, and one sets 

5- = U 5-, 5— = n sy, (39) 

a;gR a;GR 

where 5“°° is independent of p. 

Consider now a manifold M endowed with a connection. We denote by 
Ly{M) the space of operators on C°°{M) which locally are pseudodifferential 
operators with symbols (as dehned in ()18p l in Sy. Given a linear function 
l{v,x) as defined in section [3l let V’ £ C°°{M x M) such that it is 1 on a 
neighbourhood of the diagonal and such that dxl{v,x) y 0 for il!{xo,x) / 0 
and 0 / u G T* 
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Then, for any operator A G Lp, define 


cTAiv) 




x=7r 


(v) 


(40) 


It can be proved [3] that such a function belongs to S^, and that different 
choices of functions V’ and I lead to the same function a a modulo on element 
of S~°°. Therefore, the symbol a a of the operator A is defined as the 
corresponding equivalence class in Sp /S~°°. 


8 Symbols in quantum field theory 

On the side of physical applications, let us here reconsider the photon prop¬ 
agator in the Euclidean version of quantum electrodynamics. In modern 
language, the functional integral tells us that the photon propagator is ob¬ 
tained by first evaluating the gauge-field operator /r, i/ = 0,..., 3, re¬ 
sulting from the particular choice of gauge-averaging functional, then taking 
its symbol cr(P^i/) and inverting such a symbol to find cr~^{Pp,^) = for 
which (tS = Scr = I. The photon propagator reads eventually [3] 

(x, y) = (27r)-^ ^ d^k ( 41 ) 

for some contour The gauge-field Lagrangian turns out to be 

^ (42) 

where 

p^ = \AAd^,A’^-d’^Ap) + ^Apd’'A,, (43) 

and 

Note that here we denoted by the standard partial derivative with respect 
to the fi component, that is Also, a G M\{0}, gfj,^ = diag(l, 1,1,1) 

and □ = g^'^dfj,dy. Of course, the term pp only contributes to a total diver¬ 
gence and hence does not affect the photon propagator, while the parameter 
a can be set equal to 1 (Feynman choice) so that calculations are simplified. 
Thus, we can eventually obtain the gauge-field operator 

Ppy{a = 1) = -gpu'A. (45) 

Its symbol, which results from Fourier analysis of the □ operator, reads as 

cF{Ppu{a = 1)) = k^gpu, (46) 
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and hence the Euclidean photon propagator is (cf. eq. dH) 

(47) 

where the points x and y refer to the indices y and v, respectively. Note that 
integration along the real axis for ko,ki,k 2 , k^ avoids poles of the integrand, 
which are located at the complex points for which k'^ = ~ 

Strictly speaking, the gauge parameter a in (I4H]1 and (|i^ is the bare 
value as of a before renormalization, and we should express the bare symbol 
of the gauge-field operator in QED in the form 

= k^giiv + ( —-1) = CTfiuik). (48) 

\ J 

Its inverse is a combination of and k^k^ with coefficients A and B, 
respectively, determined from the condition 




which implies 


A = — 5 = 

fe2’ ' 


At this stage, the bare photon propagator takes the form 

d^k 


ix,y) = 


i2nr 




ik-{x-y) 


(49) 

(50) 


(51) 


9 Hawking radiation 


Following Ref. [5] we now consider a completely different setting, i.e. the 
spacetime of a collapsing star, for which, outside the horizon, the wave 
operator appears as it would in a Schwarzschild spacetime for a suitable 
choice of coordinates. Let r be such a radial coordinate with associated 
hypersurface S having equation r = 0, and let t be Killing time, so that 
k = -§1 IS the timelike Killing vector field outside the Killing horizon B^- 
The wave operator reads as 

_ 9^ 1 / 2 m\ d f 2 f d \ 1 / 2m 

dP r2 r y (9r \ \ r ) dr) ^ r 



where the operator L is independent of f or r and takes the same form as in 
Minkowski spacetime. The radial part of the wave operator is 


ot^ \ r J or \ \ r J Or J 


(53) 
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By letting 7 ^ (here d is for detector, since a calculation along 7 ^ yields the 
spectrum as measured by an observer whose worldline is jd) be an integral 
curve of ^ and defining t = Kt, the radial part of the wave operator can be 
decomposed in the form 


n„ = (..7 - Ah = ^ (1 - iL (1 _ a ) , (54) 

where the dot denotes 

Since the operator y/Aji is on firm ground, the operator D/j can be 
decomposed as 



K'jd + 



In the following, we will investigate the spectral properties of the operator 
K'jd — y/Aji, corresponding to the outgoing part of the radiation. Let us start 
by considering the first-order pseudodifferential operator (our conventions 
for numerical factors follow here our Ref. m) 

+ (55) 


obtained by taking the Fourier transform with respect to the variable r. 

In order to study the spectrum of \/A^, we consider the eigenvalue 
equation Aru + )^u = 0, where u depends only on r. Upon denoting by a 
prime the differentiation with respect to r, this reads as 


2m\ 


r J 


( 1 _ zl: ) ^') + ( 1 - — ) 


2m\ 


-1 


r J 


This ordinary differential equation implies that u{2m) 
finite. At large r, eq. (fSUI) reduces to 


u = 0. (56) 

= 0 while u'{2m) is 


{r^u'y + XW = 0 . 


(57) 


The general solution of eq. dSTD that is bounded at infinity is [5] 


, , sin(Ar) cos(Ar) 
n(r) = a -^^ -b b -^ 


(58) 


where the parameters a and b are constant. Prom the condition u{2m) = 0 
one hnds 

-= tan(2mA), (59) 

a 


and [5] 


An = A-h n—, ne Z 

2m 2m 


(60) 
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is the large-r limit of the spectrum of A/j, having set 6 = arctan(—6/a). 
Since the action of the wave operator on smooth functions should be smooth 
on the hypersurface S, one also requires boundedness of n as r —)• O'*". This 
implies in turn that b = 0 = 0 and Xn reduces to n upon rescaling the radial 
variable r. 

The spectral (/-function for the operator P defined in (|55l] is therefore 
expressed, in the large-r limit, by the asymptotic expansion 

CXD 

n=l 


where the tilde is used to denote removal of the degeneracy of the vanishing 
eigenvalue. By relabelling the lower limit of summation, we can re-express 
this spectral ^-function in terms of the Hurwitz (/-function, i.e. [5] 

OO 

Cp('S) ~ ^ + 1 + 1 —= Cp (s, 1 + 1—• (62) 

1=0 


If uj is the variable dual to ^ in the framework of Fourier transform, one 
can write that the Fourier transform with respect to ^ of Cp{s) is approxi¬ 
mated by [6] 


^{Cp{s)) ~ 




S o_1 


F(s)(e2W« - 1)' 


(63) 


If one takes the limit as s —)• 1 and recalls, from Ref. [7], that k = ^ in the 
spacetime of a collapsing star, one finds 


/ XN Sttui , ,, 

The spectral density p is now given by 

~ Stt 

p = kP{Cp{s)) ~ , (65) 

which is the famous result of Hawking in Ref. [8] . 


10 Concluding remarks 

Pseudodifferential and Fourier-Maslov integral operators [10] play a key 
role in the modern theory of elliptic and hyperbolic equations on mani¬ 
folds, respectively, and the physical applications form an equally rich family, 
ranging from the Cauchy problem of classical field theory m to the Green 
functions of quantum held theory and black-hole physics, as we have shown. 

Here we would like to add that, in the sixties, DeWitt discovered that 
the advanced and retarded Green functions of the wave operator on metric 
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perturbations in the de Bonder gauge make it possible to define classical 
Poisson brackets on the space of functionals that are invariant under the 
action of the full diffeomorphism group of spacetime. He therefore tried 
to exploit this property to dehne invariant commutators for the quantized 
gravitational field [12], but the operator counterpart of the classical Poisson 
brackets turned out to be a hard task. On the other hand we know from 
section m that, rather than inverting exactly a partial differential opera¬ 
tor, it is more convenient to build a parametrix. This makes it possible to 
solve inhomogeneous equations with the desired accuracy. Interestingly, it 
remains to be seen whether such a construction might be exploited in canon¬ 
ical quantum gravity, provided one understands what is the counterpart of 
classical smoothing operators in the quantization procedure. 
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A Symbol map for partial differential operators; 
the space of pseudodifferential operators on 

The definitions of Section 1 are best suited to deal with the analysis of 
inhomogeneous partial differential equations Pu = / and use a nomenclature 
very close to the one appropriate for Fourier-Maslov integral operators for 
hyperbolic equations. However, we should also recall the basic properties 
summarized below m- 

A linear partial differential operator P of order d on is a polynomial 
expression 

P{x,D) = ^ a„(x)B“, a„ E (66) 

(y:\a\<d 

where, for the multi-index a = (ai,..., a„), the modulus |a| and the deriva¬ 
tive operator are defined by 

/ o \ «i / f) 

H^o. + ... + a„, i)J^H)H(_) ...(_) . (67) 

The symbol of P is then defined by 

a{P) = a{x,^)= y] aa(a^)C", (68) 

Q::|Q:|<d 
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and is a polynomial of order d in the dual variable where (x,^) defines a 
point of the cotangent bundle of M”. The leading symbol is the highest order 
part of (T(a:, i.e. 


crL(P) = ^ aaix)C, (69) 

a: |Q:|=d 

and the action of P can be re-expressed in integral form as 

Pf{x) = {2Tr)-^ [ e‘^'V(x,0/(OdC- (70) 

In general, one can consider the set S'^ of all symbols a{x,^) such that 

(i) a is smooth in (x,^) with compact x support. 

(ii) For all multi-indices {a,j3), there exist constants Ca,p for which 


D^Dlaix,0 <C„,^(1 + |^|) 


d-l 


(71) 


For a E S'^, one dehnes the associated operator (5 being the Schwartz 
space of smooth complex-valued functions with fast decrease) P : 5 —)• 
C^(M”) as in (A5), i.e. 

Pf{x) = (27r)"t f e‘^'V(x,C)/(0d^ = (2vr)"” [ P^'^^~y'^a{x,C)f{y)dydC 

7R" JR" 

(72) 

The space of such operators is the set of pseudo-differential operators of 
order d In- 
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